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The possibility of ferroelectricity in rubidium niobate 
and rubidium tantalate with the perovskite structure was 
discussed by Smolenskii and Kozhevnikova^ and then by 
Megaw^ in the early 1950s. In Ref. Q1 the authors re¬ 
ferred to unpublished data by V.G. Prokhvatilov who 
detected the tetragonal RbTa0 3 phase with a = 3.92 A, 
c = 4.51 A exhibiting a phase transition near 520 K; 
in Ref. @ these data were simply cited. However, fur¬ 
ther studies have shown that, unlike lithium, sodium, and 
potassium niobates, RbNb0 3 and RbTa0 3 crystallize in 
individual crystal structures with triclinic PI symmetry 
for RbNb0 3 and monoclinic C2/m one for RbTaO^— 
when prepared at atmospheric pressure. To obtain these 
materials with the perovskite structure, they should be 
prepared at high pressures (65-90 kbar)j£ Due to the dif¬ 
ficulties in synthesis of RbNb0 3 and RbTa0 3 with the 
perovskite structure, the properties of these crystals have 
been studied very little. 

The phase diagrams of Rb20-Nb20,5 and Rb20 
Ta20.5 systems were studied in Ref. 0 and0. RbNb0 3 is 
formed by the peritectic reaction and decomposes above 
964° Cj£ RbTa0 3 decomposes above 600° C probably due 
to the peritectic reaction too£ Rubidium-containing fer¬ 
roelectric materials in the BaNb 206 -NaNb 0 3 -RbNb 0 3 
system with the tungsten-bronze structure have high 
electro-optical properties that substantially exceed those 
of lithium niobate^'iS In Ref. the possibility of us¬ 
ing rubidium niobate and rubidium tantalate for pho¬ 
toelectrochemical decomposition of water was discussed. 
In Ref. [HI, it was proposed to use the delamination of 
RbTa0 3 structure to produce porous Ta0 3 nanomem¬ 
branes with pore sizes of 1.3 x 0.6 and 1.1 x 1.1 A, which 
can be used for selective filtration of lithium ions. 

The lack of knowledge on the properties of compounds 
under consideration appears, in particular, in contradic¬ 
tory data on the ferroelectric properties of RbTa0 3 . For 
example, the existence of the phase transition at 520 K 
in the tetragonal phase was reported in Ref. T, whereas 
the data of Ref. [1 showed that RbTa0 3 prepared at high 
pressure has the cubic perovskite (or close to it) struc¬ 


ture. At 300 K, the structure of RbNb0 3 is similar to 
that of the orthorhombic BaTi0 3 , and the data of differ¬ 
ential thermal analysis indicate phase transitions in it at 
15, 155, and 300°C£ 

In this work, the equilibrium structures of RbNb0 3 
and RbTa0 3 were determined from first-principles cal¬ 
culations, and spontaneous polarization, dielectric con¬ 
stant, piezoelectric and elastic moduli, nonlinear optical 
and electro-optical properties as well as the energy band 
gaps in the LDA and GW approximations were calcu¬ 
lated for these crystals. 

The first-principles calculations were performed within 
the density functional theory using the ABINIT soft¬ 
ware^ The exchange-correlation interaction was de¬ 
scribed in the local density approximation (LDA). The 
optimized norm-conserving pseudopotentials for Nb, Ta, 
and O atoms used in these calculations were taken from 
Ref. HD- The non-relativistic pseudopotential for the 
Rb atom (electronic configuration 4s 2 4jp 6 5s°) was con¬ 
structed according to the scheme of Ref. (l5| using the 
OPIUM program-^ with the following parameters: r s = 
1.68, r p = 1.72, r d = 1.68, q s = 7.07, q p = 7.27, 
qd — 7.07, r m in = 0.01, r ma x — 1.52, and Vioc = 1-58 a.u. 
(for notations, see Ref. G3). The testing of the Rb pseu¬ 
dopotential on the PI phase of RbNb0 3 and the C2/m 
phase of RbTa0 3 , which are stable at atmospheric pres¬ 
sure, showed its sufficiently high quality: the calculated 
lattice parameters and atomic coordinates in these phases 
(see Tables Q] and [Til are in good agreement with the 
experimental dataj^ small underestimates of the calcu¬ 
lated lattice parameters are characteristic of the LDA 
approximation used in this work. 

The lattice parameters and equilibrium atomic posi¬ 
tions in the unit cells were determined from the condition 
when the residual forces acting on the atoms were be¬ 
low 5T0 -6 Ha/Bohr (0.25 rneV/A) in the self-consistent 
calculation of the total energy with an accuracy better 
than 10~ 10 Ha. The maximum energy of plane waves 
was 30 Ha for RbNb0 3 and 40 Ha for RbTa0 3 . Inte¬ 
gration over the Brillouin zone was performed using a 
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TABLE I. Calculated lattice parameters and atomic coordi¬ 
nates in RbNbOs structures. 


Atom 

Position x 

y 

z 



Phase PI 




a = 

5.0816, 6 = 8.3047, c 

= 8.7916 A, 



a = 114.0625, 0 = 93.3891 

7 = 95.1160° 


Rb2 

la 

+ 0.00000 

+ 0.00000 

+ 0.00000 

Rbl 

ib 

+ 0.00000 

+ 0.00000 

+0.50000 

Rb3 

2 i 

+0.41251 

+0.70257 

+0.09488 

Nbl 

2 i 

+0.49674 

+0.28138 

+0.35602 

Nb2 

2 i 

+0.02746 

+0.51037 

+0.30988 

Ol 

2 i 

+0.10125 

+0.39078 

+0.82160 

02 

2 i 

+0.23664 

+0.42747 

+0.50994 

03 

2 i 

+0.28650 

+0.71922 

+0.45293 

04 

2 i 

+0.29777 

+0.37205 

+0.19910 

05 

2 i 

+0.33951 

+0.05579 

+0.27137 

06 

2 i 

+0.78420 

+0.27571 

+0.22313 



Phase Pm3m 




a = 4.0291 A 


Rb 

la 

+ 0.00000 

+ 0.00000 

+ 0.00000 

Nb 

lb 

+0.50000 

+0.50000 

+0.50000 

O 

3c 

+ 0.00000 

+0.50000 

+0.50000 



Phase P4mm 




a = 4.0037, c = 4.1592 A 


Rb 

la 

+ 0.00000 

+ 0.00000 

-0.00336 

Nb 

lb 

+0.50000 

+0.50000 

+0.51818 

Ol 

2c 

+0.50000 

+ 0.00000 

+0.47446 

02 

lb 

+0.50000 

+0.50000 

-0.03654 



Phase Amm2 



a = 

3.9928, 6 = 5.7742, c 

= 5.7960 A 


Rb 

2 a 

+ 0.00000 

+ 0.00000 

-0.00341 

Nb 

26 

+0.50000 

+ 0.00000 

+0.51447 

Ol 

4e 

+0.50000 

+0.25496 

+0.22842 

02 

2 a 

+ 0.00000 

+ 0.00000 

+0.47735 



Phase R3m 





a = 4.0571 A, a = 89.8945° 


Rb 

la 

-0.00308 

-0.00308 

-0.00308 

Nb 

la 

+0.51193 

+0.51193 

+0.51193 

O 

36 

-0.02115 

+0.48415 

+0.48415 


8 x 8 x 8 Monkhorst- Pack mesh. The spontaneous polar¬ 
ization in ferroelectric phases was calculated by the Berry 
phase method. The phonon spectra, dielectric constants, 
piezoelectric and elastic moduli were calculated within 
the density-functional perturbation theory similarly to 
Ref.[H Nonlinear optical and electro-optical properties 
were calculated using the technique described in Ref. HH. 
All physical properties presented in this work were cal¬ 
culated for the theoretical lattice parameter. 

The phonon spectrum of RbNbC >3 in the cubic Pm3m 
phase is shown in Fig. 1. This spectrum contains a band 


TABLE II. Calculated lattice parameters and atomic coordi¬ 
nates in RbTa 03 structures. 


Atom 

Position 

X 

y 

Z 



Phase C2/m 




a = 

b = 6.3396, c = 

8.0171, 



a = 86.1031, 0 = 93.8969, 7 = 96.8997° 


Rbl 

4l 

+0.16000 

-0.16000 

+0.73758 

Rb2 

4 g 

+0.26494 

+0.26494 

+ 0.00000 

Tal 

4 h 

+0.31104 

+0.31104 

+0.50000 

Ta2 

4 i 

+0.23924 

-0.23924 

+0.30214 

Ol 

8 j 

+0.27303 

+0.04639 

+0.39643 

02 

8 j 

+0.54749 

-0.21209 

+0.28704 

03 

4 i 

+0.16752 

-0.16752 

+0.08653 

04 

4 i 

+0.37705 

-0.37705 

+0.55235 



Phase Pm3m 




a = 3.9846 A 



Rb 

la 

+ 0.00000 

+ 0.00000 

+ 0.00000 

Ta 

16 

+0.50000 

+0.50000 

+0.50000 

O 

3c 

+ 0.00000 

+0.50000 

+0.50000 
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FIG. 1. Phonon spectrum of RbNbOs in the cubic Pm3m 
phase. Labels near curves indicate the symmetry of unstable 
modes. 


of unstable modes characteristic of ferroelectric chain in¬ 
stability which was first observed in KNbC> 3 ^ At the 
center of the Brillouin zone, this mode has the T 15 sym¬ 
metry, is triply degenerate, and describes the ferroelectric 
distortion of structure. The structures appearing upon 
condensation of the A '5 and M 3 modes are characterized 
by antiparallel orientation of polarization in neighboring 
...-O-Nb-O-... chains. 

The energies of all RbNbOs phases formed upon con¬ 
densation of the above unstable modes are given in Ta¬ 
ble IIIII Among these phases, the R3m phase has the 
lowest energy. The phonon spectrum calculations for 
the i?3m phase show that the frequencies of all optical 
phonons at the center of the Brillouin zone and at high- 
symmetry points at its boundary are positive; the deter- 













































3 


TABLE III. Relative energies of low-symmetry RbNbOs 
phases formed from the cubic perovskite phase upon conden¬ 
sation of unstable phonons, phases with 6 H, AH, 9 R, and 2H 
structures, and the Pi phase prepared at atmospheric pres¬ 
sure (the most stable phase energy is in boldface). 

Phase 

Unstable mode Energy, meV 

Pm3m 

— 

0 

PA/nmm 

Mi 

-31.3 

Pmma 

x 5 

-34.8 

Cmcm 

X 5 

-38.7 

PAmm 

r 15 

-46.5 

Amm2 

r 15 

-57.0 

R3m 

r 15 

- 58.6 

PI 

— 

+27.4 

PQs/mmc (6 H) 

— 

+ 121.4 

P&z/mmc (AH) 

— 

+334.0 

R3m (9 R) 

— 

+568.1 

PQs/mmc (2 H) 

— 

+1752 


TABLE IV. Relative energies of low-symmetry RbTaOs 
phases formed from the cubic perovskite phase upon conden¬ 
sation of unstable phonons, phases with 6 H, AH, 9 R, and 
2 H structures, and the C2/m phase prepared at atmospheric 
pressure (the most stable phase energy is in boldface). 

Phase 

Unstable mode 

Energy, meV 

Pm3m 

— 

0 

PAmm 

r 15 

-1.80 

Amm2 

r 15 

-1.87 

R3m 

r 15 

- 1.90 

C2/m 

— 

+38.5 

P& 3 /mmc (6 H) 

— 

+113.3 

PQz/mmc (AH) 

— 

+352.2 

R3m (9 R) 

— 

+589.0 

PQz/mmc (2H) 

— 

+1839 


minant and all leading principal minors constructed of 
elastic moduli tensor components are also positive. This 
means that the R3m phase is the ground-state struc¬ 
ture of RbNbC> 3 . The calculated lattice parameters and 
atomic coordinates in this phase are given in Table |I] As 
the same sequence of phases as in BaTiC >3 is supposed in 
rubidium niobate with the perovskite structure^ the lat¬ 
tice parameters and atomic coordinates in two other fer¬ 
roelectric phases are also given in this table. The lattice 
parameters calculated for the orthorhombic RbNbC >3 are 
in good agreement with the experimental data obtained 
at 300 K (a = 3.9965, b = 5.8360, and c = 5.8698 A, 
Ref.@). 

In RbTaC> 3 , the frequency of unstable Tis phonon in 
the phonon spectrum (Fig. 2) and the energy gain re¬ 
sulting from the transition to ferroelectric phases (Ta¬ 
ble [TV]) are rather low; so it is necessary to additionally 
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FIG. 2. Phonon spectrum of RbTaC >3 in the cubic PmSm 
phase. Labels near curves indicate the symmetry of unstable 
modes. 


test the stability of the ferroelectric distortion with re¬ 
spect to zero-point lattice vibrations. For this purpose, 
we used the technique proposed in Ref. HI The energy 
gain resulting from the transition from the Pm3m phase 
to the R3m phase is Eq = 1.90 meV and the unstable 
phonon frequency at the T point in the Pm3m phase is 
v = 84 cm -1 . As the energy ratio his/E 0 = 5.51 exceeds 
the critical value of 2.419 obtained in Ref. HI the energy 
of the lowest vibrational state in a two-well potential ap¬ 
pears above the upper point of the energy barrier sepa¬ 
rating the potential wells, and the ferroelectric ordering 
is suppressed by zero-point vibrations. Therefore, the 
only stable phase of RbTaOs with the perovskite struc¬ 
ture is the cubic phase. The calculated lattice parameter 
of this phase is given in Table El its value is in satisfac¬ 
tory agreement with the experimental data (a = 4.035 A, 
Ref. [H). 

It is known that the formation of phases with hexag¬ 
onal BaNiC >3 (polytype 2 H), hexagonal BaMnOs (poly¬ 
type 4 H), hexagonal BaTiC >3 (polytype 6 H), and rhom- 
bohedral BaRu 03 (polytype 9 R) structures is character¬ 
istic of AB O 3 perovskites with the tolerance factor t > 1, 
and the studied compounds belong to this class. Our 
calculations showed that the energies of these phases for 
both rubidium compounds are appreciably higher than 
the cubic phase energy f Tables Hill and ITVl) . These results 
explain why it was impossible to observe the transition 
of RbNbOs to the hexagonal structure^ upon heating, by 
analogy to that occurring in BaTiOs- The high energies 
of these phases, in particular, the 2H phase, are probably 
caused by larger sizes and strong electrostatic repulsion of 
Nb 5+ ions which occupy face-sharing octahedra in these 
structures. 

We consider now some properties of ferroelectric 
RbNbC> 3 . The calculated polarization in RbNbC >3 is 
0.46 C/m 2 in the PAmm phase and 0.50 C/m 2 in the 
Amm2 and R3m phases; these values slightly exceed 
the calculated polarization in the same phases of KNbOs 
(0.37, 0.42, and 0.42 C/m 2 , respectively). The static di- 
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TABLE V. Nonzero components of the piezoelectric tensor a v 
(C/m 2 ) and tensors of the second-order nonlinear optical sus¬ 
ceptibility di„ and the linear electro-optic effect r (pm/V) 
in rhombohedral phases of RbNbOs, KNbC> 3 , LiNbOa, and 
BaTiOs. 


Coefficient 

RbNbOs 

KNbOs 

LiNbOs 

BaTiOs 

en 

-3.0 

-4.2 

-2.4 

-4.0 

ei5 

+4.8 

+6.8 

+3.5 

+7.3 

e3i 

+2.4 

+2.3 

+0.1 

+3.5 

e33 

+2.9 

+3.1 

+1.1 

+5.1 

dn 

+ 12.7 

+11.9 

+2.3 

+4.4 

dis 

-23.6 

-21.9 

-11.5 

-16.1 

d3i 

-23.6 

-21.9 

-11.5 

-16.1 

d33 

-29.4 

-27.3 

-37.4 

-31.1 

rn 

-12.8 

-17.7 

-5.6 

-13.7 

ns 

+27.6 

+39.2 

+17.1 

+43.3 

+31 

+18.0 

+23.9 

+10.1 

+25.3 

r33 

+30.1 

+40.6 

+27.3 

+48.9 


electric tensor in the R3m phase is characterized by two 
eigenvalues: ejj = 21.1 and £*}_ = 35.8; the optical dielec¬ 
tric tensor eigenvalues are £jj° = 5.31 and £^° = 5.91. 
In the cubic phase, the elastic moduli are C\\ = 412, 
C12 = 84, and C 44 = 102 GPa; the bulk modulus is 
B = 193.5 GPa. The nonzero components of the tensors 
of piezoelectric effect e,;„, second-order nonlinear optical 
susceptibility di V . and linear electro-optic (Pockels) effect 
Ti V in the R3m phase of rubidium niobate are compared 
with the corresponding properties of other rhombohedral 
ferroelectrics in Table El We can see that the piezoelec¬ 
tric moduli in rhombohedral RbNbOa (as well as in its 
other polar phases) are slightly lower than in KNbOa. 
The nonlinear optical coefficients in RbNbOs exceed the 
corresponding values in KNbC> 3 , although the ^33 value 
in rubidium niobate is slightly lower than in lithium nio¬ 
bate. As for the electro-optical properties, in rhombo¬ 
hedral RbNbC >3 they are slightly lower than in KNbC> 3 , 
but are notably superior to those of lithium niobate. In 
the orthorhombic phase (stable at 300 K), nonlinear opti¬ 
cal properties of RbNbC >3 are comparable to those of the 
same phase of potassium niobate: for example, the ^33 
modulus is —30.8 pm/V in RbNbC >3 and —30.4 pm/V in 
KNb0 3 . 

In cubic RbTaOa, the optical dielectric constant is 
£oo = 5.58. The static dielectric constant can be es¬ 
timated only in the rhombohedral phase as ~140. The 
elastic moduli in cubic rubidium tantalate are CVi = 466, 
Ci 2 = 91.5, and C44 = 120 GPa; B = 216 GPa. The 
piezoelectric moduli, second-order nonlinear optical sus¬ 
ceptibility, and electro-optical coefficients in the cubic 
phase are zero. 

An unexpected result of our calculations is that in 
both studied compounds the PI and C 2 /m phases which 
can be prepared at atmospheric pressure are metastable. 


This result is probably caused by an effective lattice con¬ 
traction which always exists in the LDA calculations. 
The fact that the specific volume of the Pm3m phase 
is noticeably smaller than that of Pi, C2/m, P&z/mmc, 
and R3m phases suggests that under pressure the cubic 
perovskite phase will be the most stable one. To esti¬ 
mate the maximum value of the actual effective pressure, 
the lattice parameters and atomic positions in the C2/m 
structure of rubidium tantalate were calculated for differ¬ 
ent pressures and it was shown that the unit cell volume 
equal to the experimental one at 300 K can be obtained 
at an isotropic pressure of —24.7 kbar. At this pressure, 
the enthalpy of the C2 /to phase becomes lower than that 
of Pm3m by ~230 meV, i.e., becomes consistent with 
the experimental data. At the above-mentioned nega¬ 
tive pressure, the ratio hv / Eq determining the stability 
of the ferroelectric phase in RbTaC >3 becomes equal to 
1.90, i.e., it is slightly lower than the critical value of 
2.419. However, if we take into account that the above 
negative pressure is obviously overestimated because it 
includes the thermal expansion effect, we can suppose 
that, even taking into account the systematic error in 
the LDA lattice parameter determination, rubidium tan¬ 
talate will remain cubic up to the lowest temperatures. 

The conclusion that RbTaOs is an incipient ferroelec¬ 
tric in which the ferroelectric ordering is suppressed by 
zero-point vibrations agrees with the data of Ref. but 
contradicts the data of Ref. |H in which the phase transi¬ 
tion near 520 K was reported. We suppose that tantalum- 
enriched phases (in particular, with the tungsten-bronze 
structure^) could be formed in rubidium tantalate sam¬ 
ples discussed in Ref. [l] because of the low temperature 
of the peritectic reaction, and this could result in the 
observed anomaly. 

In Ref. the possibility of using various oxides 
with the perovskite structure, in particular RbNbOs and 
RbTaC> 3 , for development of photoelectrochemical solar 
cells was discussed. We calculated the band gap E g in 
these compounds both in the LDA approximation and 
in the GW approximation that takes into account many- 
body effects (the technique of the latter calculations was 
analogous to that used in Refs. l2ll - [23l ). In the LDA ap¬ 
proximation, E^ da = 1.275 eV in cubic RbNbOs when 
the spin-orbit coupling is neglected; in P4 toto, Amro2, 
and R3m phases, E g DA is 1.314, 1.869, and 2.137 eV, 
respectively. In cubic RbTaOs, E g DA = 2.175 eV when 
the spin-orbit coupling is neglected. The valence band 
extrema in the cubic phase of both compounds are at the 
R point of the Brillouin zone, whereas the conduction 
band extrema are at the T point. The calculations using 
the technique of Ref. [23] yield the spin-orbit splitting of 
the conduction band edge A 50 = 0.111 eV for RbNbOs 
and A so = 0.400 eV for RbTaOs; the spin-orbit splitting 
of the valence band edge is absent. After correction for 
the conduction band edge shift (Ago/3), the LDA values 
of E g that take into account the spin-orbit coupling are 
1.238, 1.277, 1.832, 2.100, and 2.042 eV for four RbNb0 3 
phases and for cubic RbTaOs, respectively. 
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In the GW approximation, the band gap when the 
spin-orbit coupling is neglected is E^ w = 2.403, 2.616, 
3.291, and 3.609 eV, respectively in cubic, tetragonal, 
orthorhombic, and rhombohedral RbNbOa and 3.302 eV 
in cubic RbTaOa- If the spin-orbit coupling is taken 
into account, these values decrease to 2.366, 2.579, 3.254, 
3.572, and 3.169 eV, respectively. The obtained values 
of E a are appreciably smaller than those calculated in 
Ref. [ll] for cubic phases (3.4 eV for RbNbC >3 and 4.3 eV 
for RbTa0 3 ). 

Some authors who studied rubidium niobate and ru¬ 
bidium tantalate have noticed their sensitivity to humid¬ 
ity. Evidently, this can be a serious obstacle for practi¬ 
cal applications of these materials. However, we would 
like to note that this property is inherent to phases pre¬ 
pared at atmospheric pressure and having loose struc¬ 
tures whose specific volume is 26-28% larger than that 
of the perovskite phase. In Ref. [§, it was suggested that 


the effect is due to intercalation of water molecules into 
the loose structures, rather than to hydrolysis of these 
compounds. The possibility of preparing RbTaOa by hy¬ 
drothermal synthesis 24 ' and the low rate of RbTaC >3 ion 
exchange in HC1 during its delaminationi 2 support this 
idea. This suggests that the considered compounds with 
the perovskite structure can be quite stable to humidity. 

In summary, the present calculations of RbNbC >3 and 
RbTaC >3 properties and their comparison with the prop¬ 
erties of other ferroelectrics show that rubidium niobate 
is an interesting ferroelectric material with high nonlin¬ 
ear optical and electro-optical properties, and rubidium 
tantalate is an incipient ferroelectric. 

The calculations presented in this work were performed 
on the laboratory computer cluster (16 cores). 

This work was supported by the Russian Foundation 
for Basic Research, project no. 13-02-00724. 
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